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The intra-atomic magnetic dipole moment - frequently called (T^) term - plays an important role in the 
determination of spin magnetic moments by x-ray absorption spectroscopy for systems with nonspherical 
spin density distributions. In this work, we present the dipole moment as a sensitive monitor to changes in 
the electronic structure in the vicinity of a phase transiton. In particular, we studied the dipole moment at 
the Fe^"^ and Fe^^ sites of magnetite as an indicator for the Verwey transition by a combination of x-ray 
magnetic circular dichroism and density functional theory. Our experimental results prove that there exists 
a local change in the electronic structure at temperatures above the Verwey transition correlated to the 
known spin reorientation. Furthermore, it is shown that measurement of the dipole moment is a powerful 
tool to observe this transition in small magnetite nanoparticles for which it is usually screened by blocking 
effects in classical magnetometry. 

X-ray absorption near edge spectroscopy (XANES) and its associated magnetic circular dichroism (XMCD) 
is an experimental technique that allows for the determinination of magnetic properties element- specif- 
ically. Today, it is widely spread as one main advantage of XMCD that the orbital moment and the spin 
moment can be determined independently from each other using sum rules ^'^ which have been experimentally 
confirmed for the 3d transition metals^. The intra-atomic magnetic dipole moment of the 3d spin-density 
distribution is one term in the XMCD sum rule for the spin moment. More precisely, with the spin sum rule 
the effective spin moment {lds,eff), which is the sum of the intrinsic spin moment (—2 (Sz) iib) and the dipole 
moment (7 (T^) plb) of the spin density distribution, can be determined from the integrals of the XMCD spectrum 
and the isotropic absorption spectrum^. Since the sum rule contains two terms, the contributions of the intrinsic 
spin moment and the magnetic dipole moment cannot be readily determined. The latter is known to be negligibly 
small for cubic symmetry only^. Limitations of the spin sum rule and the importance of the magnetic dipole term 
in noncubic environments such as surfaces and interfaces were studied with electronic band structure calculations 
for iron, cobalt and nickel^. The seperate determination of spin and magnetic dipole moments with an angle 
averaging spin sum rule was demonstrated theoretically^ and experimentally^ for a highly anisotropic 3d trans- 
ition metal system where the charge and spin densities were not spherically symmetric. The magnetic dipole 
moments of single cobalt atoms and nanoparticles (NPs) were quantified by combining the effective spin moment 
determined from measured XMCD spectra using the spin sum rule with the intrinsic spin moment calculated in 
the local spin density approximation^. According to ah initio studies, the (T^) term, in general, cannot be neglected 
in the spin sum rule for systems with low dimensionalities like monolayers and monatomic wires^ and clusters^. 
Anisotropics of the spin density distributions causing large magnetic dipole moments were also investigated in 
molecules on surfaces like Cu-phthalocyanine on Ag(001)^° and Fe-octaethylporphyrin on Cu(001)^\ 

For many examples, the measured contribution of the magnetic dipole moment can be regarded as a dis- 
advantage of the XMCD impeding the extraction of the spin moment in a straight-forward way. In this work, we 
employ it as a monitor for changes in the electronic structure as it is caused by phase transitions. In order to show 
the potential of the XMCD technique, we chose magnetite (Fe304) as an example because in magnetite the average 
Fe dipole moment in the low-temperature phase according to our electronic structure calculations is sufficiently 
large to be reliably quantified with XMCD. In addition, for magnetite XMCD is also site- selective giving the 
possibility to distinguish between Fe ions at different lattice sites. 

Magnetite crystallizes in the inverse spinel structure in a cubic (isometric) crystal system, with the oxide anions 
arranged in a cubic close-packed lattice while the Fe cations are located on octahedral and tetrahedral lattice sites 
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in between. In particular, Fe^^ occupies the tetrahedral sites (also 
denoted A sites in the literature) and half of the octahedral sites (B 
sites), whereas the Fe^^ ions are located at the remaining octahedral 
sites. For a schematic diagram of the monoclinic P2/c unit cell of 
magnetite containing 56 atoms we refer to Fig. 1 in Ref. [17]. The 
strongest interaction is the superexchange between Fe^^ on octahed- 
ral and Fe^^ on tetrahedral sites which is negative resulting in an 
antiparallel orientation of spins of different sublattices. 

At low temperatures, magnetite undergoes a phase transition that 
was discovered in 1939 by Verwey as a rapid change of the electronic 
conductivity at a certain temperature and was named after him as the 
Verwey transition (VT). In the original article^^ it was described as "a 
sudden increase of the resistance at approximately 1 17 K by a factor 
of the order 100". The resistance of a magnetite sample was com- 
pared with the one of a sample which was further oxidized and 
described as "samples containing an excess of Fe203 show a much 
smaller jump in the curve at about 120 K or even merely a change in 
the temperature coefficient"^^. Many years later it was found by heat 
capacity studies^^ that with increasing deviation 3 in Fes-s04 from 
ideal magnetite stoichiometry the VT changes from first to second 
order and that the VT temperature decreases linearly with 3. 

The refined determination of the structural changes at the VT 
based on diffraction measurements ^^'^^ initiated several theoretical 
investigations with modern methods^^"^^. The most significant struc- 
tural change connected to the VT was explained in the work of H.-T. 
Jeng et al.^^. They reported that in the low- temperature monoclinic 
unit cell, a charge and orbital ordering is related to one of the Fe^^ 
ions on an octahedral lattice site being pulled inwards the cube diag- 
onal of the sub-unit cell yielding a lower energy of the system due to 
intersite Coulomb attraction with t2g orbitals of three neighboring 
Fe^^ ions. 

With resonant x-ray scattering the changes in the structural, 
charge and orbital order were studied in dependence of the temper- 
ature across the VT using certain diffraction peaks to indicate struc- 
tural, charge or orbital order The interpretation of the diffraction 
peaks was discussed and became a controversy. 

In a structural study^^ using high-resolution synchrotron x-ray 
powder diffraction data of magnetite the charge disproportion found 
by J. P. Wright et al.^^ was confirmed between the Bl and B2 sites but 



c 




B4 

Figure 1 | Part of the magnetite unit cell in the monoclinic crystal 
structure. Oxygen atoms are drawn as small red spheres. Iron ions on 
different octahedral sites are drawn as large grey (Bl), brown (B2), green 
(B3) and blue spheres (B4). The linear three-site structure B3-B1-B3 
forms a trimeron with reduced B1-B3 bond lengths. The axes of the global 
coordinate system are also shown. 



not between the B3 and B4 sites (see Fig. 1). The picture of a wide 
distribution of different local environments around the octahedral Fe 
ions, caused by the condensation of several phonon modes, was 
favored over any bimodal charge disproportion on the octahedral 
sites^^. 

Since diffraction studies of magnetite are hampered by different 
crystallite orientations and microtwinning, the full low- temperature 
superstructure has been quite recently determined with x-ray diffrac- 
tion from an almost single-domain grain of only 40 ^m size. There it 
has also been found that localized electrons are distributed over 
linear units of three Fe ions on octahedral sites which leads to the 
observed shortening of the two Fe-Fe distances. These three- Fe-units 
were viewed as quasiparticles and named "trimerons"^^. More than 
30 years earlier the fluctuations of what nowadays is called trimerons 
were observed above the VT temperature with critical diffuse neut- 
ron scattering^^'^^ and interpreted using the Yamada model based on 
electron-phonon coupling^^ and later using a pseudospin-phonon 
theory^^. 

Below the VT temperature t2g orbital order on octahedral Fe^^ 
sites has been directly found with resonant soft x-ray diffraction 
(RSXD)^^. For this method it is important to take into account that 
there is a layer at the surface of single crystals which contributes to 
the absorption but not to the diffraction of the soft x-rays^^. By 
studying the azimuthal dependence of the diffracted intensity it 
has been shown that the distortion of the 3d orbitals towards mono- 
clinic symmetry is by far larger than that of the lattice^^ Ultrafast 
melting of the charge-orbital order leading to the formation of a 
transient phase, which had not been observed in equilibrium, was 
observed with time-resolved RSXD with a free-electron laser^^. Very 
recently it was investigated on the ps time scale how trimerons 
become mobile across the VT^°. 

In magnetite NPs the VT has already been experimentally 
observed with different methods, e.g.^ magnetization measurements 
of a non-diluted system containing strong dipolar interparticle inter- 
actions^^ and magnetoresistance measurements of tunneling junc- 
tions of stacked monolayers of magnetite NPs of 5.5 nm average 
diameter for which a VT temperature of 96 K has been determined^^. 
In these studies it has been found that the VT temperature signifi- 
cantly decreases when decreasing the size of the NPs to a small 
fraction^ ^ but that it is not very sensitive to small changes of the order 
of 10% of the particle size^^. 

However, one may note that in magnetometry, the VT can only 
hardly be distinguished from a spin reorientation transition 
that occurs in magnetite at slightly higher temperatures (Tsr ~ 
132 K^^). The spin reorientation is characterized by a vanishing mag- 
netocrystalline anisotropy constant iCi, which is negative at room 
temperature. With decreasing temperature it decreases, acquiring a 
minimum at a temperature of about 250 K^^. With further decreasing 
temperature, it increases and becomes positive at Tsr. Referring to 
the vanishing value of Ki, this temperature is also called isotropy 
point in literature. For single crystalline samples, the different trans- 
ition temperatures were measured by recording the saturation iso- 
thermal remanent magnetization along a principal axis of the single 
crystaP^. For the case of an arbitrary orientation of the magnetic field 
with respect to the principal crystallographic axes, the effects of Tsr 
and the VT on the remanence could not be clearly separated. For 
polycrystalline samples it seems to be impossible to distinguish 
between the two transitions. As pointed out by Reznicek et al.^^, there 
exists a correlation between the spin reorientation and the VT. In this 
work, the authors connected the contribution to Ki which leads to 
an anomalous increasing behavior with decreasing temperature 
between 250 K and Tsr to a local charge and orbital ordering that 
becomes long-range and stable for temperatures below the VT. 

The main purpose of the present paper is to study the magnetic 
moments, in particular the magnetic dipole moment of Fe as a func- 
tion of temperature across the phase transitions in magnetite NPs 
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with a diameter of 6 nm and a 200 nm thick polycrystalUne magnet- 
ite film which serves as a reference. For this purpose a combination of 
experimental and theoretical methods, i.e. XANES, XMCD, vibrating 
sample magnetometry (VSM) and electronic structure calculations 
based on density functional theory, has been used. As it will be 
described in detail below, a transition of the dipole moment of the 
Fe 3d spin- density distribution was experimentally found well above 
the VT and theoretically explained. 

Results 

Density functional theory. Calculations of the electronic structure 
were performed based on density functional theory on cubic 
(high temperature) and monoclinic (low temperature) structures 
of magnetite. Details of the calculations are described in the 
methods section. We have considered a monoclinic P2/c structure 
with 56 atoms per unit cell obtained by using high-resolution 
neutron and x-ray powder- diffraction data by Wright et a/.^^'^l A 
part of the monoclinic unit cell with its four inequivalent octahedral 
lattice sites, labelled B1-B4, is shown in Fig. 1. The resulting magnetic 
moments are shown in table I. The Fe sites are named as shown in 
Fig. 1 and like in Ref [17]. Our calculated values of (T^) for the cubic 
structure are negligible (—10"^ p^b) whereas for the monoclinic 
structure, they are quite large, especially for the octahedral B sites. 
It is clearly observed from the calculations that the Fe^^ ion has a 
larger contribution than Fe^^ (see Table I, where we compile the 
calculated values of 7 (T^)). At the B4 site, the contribution is 
as large as —1.44 with an anti-parallel alignment to the spin 
moment. It is only partly compensated by the second largest 
contribution of 0.72 ^5 at the Bl site with a parallel alignment to 
the spin moment. Therefore, these octahedral sites substantially 
modify the effective spin moments, as detected in the XMCD 
measurements. 

We also calculated the Fe orbital magnetic moments for the mono- 
clinic phase in the a, h and c directions as defined in Fig. 1. The 
sizeable orbital moments can be as large as 0.035 ^5 along the c 
direction for Fe ions on B4 sites and 0.031 ^5 along the h direction 
for Fe ions on Bl sites. For all other sites, the orbital moments are 
almost isotropic with absolute values of 0.015 fis to 0.019 /i^. The 
signs of the orbital moments follow the spin moments. To further 
analyze the large values of (T^), we show in Fig. 2, the m/ projected 
density of states (DOS) for different octahedral sites of the mono- 
clinic phase. The reason for the high (T^) value is the full occupancy 
of the _yi orbital in the spin down channel of the Fe^^ ion on the 
B4 site. For the Fe^^ ion, this orbital is unoccupied. Similarly, a high 
value (half of that of the B4 site, and with parallel alignment to the 
spin-moment) is observed for the Fe^^ ion on the Bl site. The large 
value also for this site is coupled to the fact that the d^^z and d-^^ 
orbitals are occupied, although only partially. Again, for the Fe^^ 
ion, both these orbitals are unoccupied, and do not contribute to (T^). 

All of our calculated results can be understood from the (T^) 
matrix elements explicitly provided in the paper by Crocombette 
et alP. For Fe^^, the spin-up channel is completely filled with almost 



Table I | Charge, spin and magnetic dipole moments for 3d orbitals 
of Fe atoms at different sites in the monoclinic unit cell. Also, effec- 
tive moments (/is,eff = ~2 (S^) iiq + V (T^) ^ib) are provided. The Fe 
sites are named as in Ref. [1 7] 
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Figure 2 | Spin-resolved m; projected density of states for the monoclinic 
structure. DOS at the tetrahedral A site along with four different 
octahedral B sites are shown. The Fe sites are named as shown in Fig. 1 and 
like in Ref. [ 17] . The specific orbitals responsible for the large contribution 
in the (T^) part are indicated by arrows. The energies in the x-axes are 
plotted with respect to the Fermi level Ep. 

no occupancy in the spin-down channel and therefore, much smaller 
(TJ contributions are observed. Note that XANES and XMCD are 
sensitive to the unoccupied DOS so that in these experiments the 
magnetic moments due to the lack of unoccupied states instead of the 
existence of occupied states is probed. Also note that even for Fe^^ in 
the monoclinic phase due to the low symmetry structure compared 
to the cubic one, we have higher contribution to (T^) compared to the 
cubic one. Obviously, there is no noticeable distinction between Fe^^ 
and Fe^^ in the cubic phase as the charge ordering is non-existent in 
the metallic phase. 

X-ray absorption spectroscopy. The XANES and XMCD spectra are 
defined as the sum and difference of two spectra obtained with 
opposite helicities of the incident circularly polarized soft x-ray 
radiation, respectively. Measured XANES spectra and corresponding 
XMCD spectra of the NPs and the film sample are shown in Fig. 3 and 
Fig. 4, respectively, at the Fe L23 absorption edges for various 
temperatures. The temperature and magnetic field histories of both 
samples were identical, i.e. both samples were field- cooled to 4 K in a 
magnetic field of 0.5 T and subsequently heated stepwise to and 
measured at the temperatures indicated by the data points in Fig. 5 
in a magnetic field of 3 T. The NPs offer significant changes of the 
XANES intensity and the XMCD intensity between 50 K and 100 K 
and the film sample between 150 K and 175 K. For both samples the 
spectra are devided into two groups: Spectra for temperatures below 
the change are shown as blue lines, and above as red lines. 

XMCD has the advantage that spin and orbital moments can be 
determined element-specifically by use of sum rules^'^. For the deter- 
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Figure 3 | Measured XANES (a) and corresponding XMCD spectra (b) of 
the nanoparticles at the 1.2,3 absorption edges. The XANES spectra are 
shown after subtracting a hnear background and normahzing at a photon 
energy in the post-edge region. 

mination of the magnetic moments via sum rules, we normaUzed the 
spectra in the usual way, i.e. first subtracted a constant background, 
then normalized at a photon energy above the absorption edge and 
finally subtracted a step like function. Please note that the spectra 
have been measured from 680 eV to 780 eV in order to be able to 
perform a reliable normalization. In Fig. 3 and Fig. 4 only the phys- 
ically relevant spectral range from 705 eV to 730 eV has been shown 
for the sake of clarity. The effective spin moment {\is,eff ^ ~2 (S^) 
+ 7 (T^) [i^ is defined as the intrinsic spin moment ( — 2 (S^) [i^ plus 
the magnetic dipole moment (7 (T^) [i^ of the spin-density distri- 
bution^ where z is defined by the magnetization direction. Its nor- 
malized values determined from the measured XMCD spectra as a 
function of temperature are shown for the NPs in Fig. 5a and for the 
film in Fig. 5b. The effective Fe spin moment of the NPs increases 
between 50 K and 100 K by about 7% relative to the value at 100 K 
and for the film between 150 K and 175 K by 16% relative to the 
value at 175 K. The magnetic moments of the NPs and the film 
sample were also measured with VSM and are shown for comparison 
as gray lines in Fig. 5a and Fig. 5b, respectively. Increases of the 
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Figure 4 | Measured XANES (a) and corresponding XMCD spectra (b) of 
the film sample at the I23 absorption edges. The XANES spectra are 
shown after subtracting a linear background and normalizing at a photon 
energy in the post-edge region. 
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Figure 5 | Normalized effective Fe spin moments as determined from the 
measured XMCD spectra (solid symbols) and the magnetic moment 
measured with VSM (thick gray line) as a function of temperature for the 
magnetite NPs (a) and the film sample (b). Errors due to the 
reproducibility of the evaluation of the XMCD spectra are within the 
symbol size. 

magnetic moments with increasing temperatures were not observed 
with VSM. This important result is due to the fact that the magnetic 
dipole moment is observable with XMCD but not with VSM as 
discussed in the following section. 

Discussion 

The electronic structure calculations show a strong change of the 
magnetic dipole moment in magnetite when the structure of the 
lattice changed from cubic to monoclinic. These changes are most 
significant for the case of Fe^^ ions on octahedral lattice sites (Bl, B4) 
which is in agreement to the experimental XMCD data. The trans- 
ition observed with XMCD (Fig. 3b and Fig. 4b) is accompanied by a 
change of the asymmetry, which is defined as the XMCD spectrum 
devided by the XANES spectrum after adequate background sub- 
traction and normalization, at the photon energies 708.2 eV and 
710 eV which mainly correspond to octahedral sites^^ The transition 
is also observed in the effective Fe spin moment (Fig. 5). However, in 
general this transition could be assigned to a change in the spin 
moment, the magnetic dipole moment or both at the same time. 
To analyse which contribution to the effective spin moment is mainly 
responsible for its significant change, the XMCD results are com- 
pared to magnetometry data obtained with VSM. The (T^) term only 
shows up in XMCD measurements because they involve an excita- 
tion from a core level to an unoccupied valence level. It does not show 
up in magnetometry or other probes that detect the macroscopic 
magnetization. For further discussion, see Ref. [39] and Ref. [40]. 
As shown in Fig. 5, there is no corresponding increase of the total 
magnetic moment in the VSM data although this method certainly is 
sufficiently sensitive to detect a change of about 7% for the NPs and 
16% for the film. Therefore the transition observed with XMCD 
cannot be due to a change in the intrinsic Fe spin moment but must 
be due to a change in the second term of the sum rule for the effective 
spin moment^, i.e. the (T^) term which describes the dipole moment 
of the Fe spin-density distribution. This conclusion is in accordance 
with our theoretical treatment of the VT. 

In Fig. 6, we compare the experimentally determined values of 
the effective spin magnetic moment of the NPs with the theoretically 
observed ones. The calculated 0 K effective moments, averaged 
over all Fe sites have been multiplied with the Bloch law 1 — 0.96 
(T/850K)^^^ to obtain magnetic moments as a function of tempera- 
ture. The theoretical moments are seen to drop at the known VT 
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temperature for the bulk, where the crystal structure changes from 
cubic to monoclinic. Again, it is the value of (T^) that changes. The 
relative change of the moments nicely agrees between theory and 
experiment in particular for the NPs although (T^) was calculated 
along one crystallographic direction (the c direction as indicated in 
Fig. 1) whereas in the measurement it was averaged over all direc- 
tions due to the random orientation of the NPs. Nevertheless, it is 
clear that a detailed measurement of (T^) shows up in an unpreced- 
ented way. This then becomes a valuable characterization tool, that 
signals not only changes in the local geometry but also the local 
electronic structure. It is likely that in general, detailed knowledge 
of (Tz) is a valuable tool for understanding local changes in the 
crystal- and electronic structure, not only for complex oxides but 
also complex actinide and rare-earth compounds. 

For 3d transition metals it has been pointed out by Stohr and 
Konig^ that spin- orbit correction terms in the perturbative calcula- 
tion of (T^) can be neglected to a good approximation, leading to the 
numerical relation (T^) + (T^) + (T^) ~ 0 and to the possibility to 
average out the magnetic dipole moment by magnetizing and mea- 
suring the sample along the three cartesian axes (a = x, y, z). In the 
case of our magnetite NPs the average magnetic dipole moment 
would vanish in each single measurement geometry because the 
NPs are randomly oriented like in a powder sample so that the 
integration of the z component = 5^(1 — 3cos^0)l2^ over the solid 
angle yields zero. However, according to Ederer et al. the effects of 
spin-orbit coupling are larger for low-dimensional systems^ and in 
particular the relation (T^) + (T^) + (T^) ~ 0 is strongly violated for 
monoatomic Fe, Co and Ni wires^\ In addition, spin-orbit coupling 
is larger for 3d oxides like magnetite than for 3d metals due to a 
stronger locaUzation of the 3d electrons. The large spin- orbit coup- 
ling becomes noticeable by the large magnetic anisotropy constant in 
the monoclinic phase of magnetite^^. Moreover, it causes that the 3d 
charge distribution is no longer independent from the spin dir- 
ection^. Therefore the effect of spin-orbit coupling in the determina- 
tion of the components {T^ can no longer be neglected, i.e. the terms 
containing the spin-flip operators do no longer vanish (see Appendix 
B in Ref. [40]). For the experimental determination of the magnetic 
dipole moment this means that at least a part of the magnetic dipole 
moment is fixed to the magnetization direction but not to the crystal 
lattice when the crystal (here each NP) is rotated relative to the 
magnetic field. In this case it is not possible to separate the intrinsic 
spin moment and the magnetic dipole moment solely by angle 
dependent XMCD measurements because the precondition for aver- 
aging out the magnetic dipole moment is no longer fulfilled. The 
present results for magnetite NPs with 6 nm diameter demonstrate 
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that the magnetic dipole moment indeed is not averaged out in the 
monoclinic low-temperature phase. 

The small average Fe orbital moments of magnetite with values 
below 0.03 ± 0.02 as determined with XMCD experiments in Ref 
[37] do not exclude sizeable local spin-orbit effects because, accord- 
ing to our calculations, the absolute values of the Fe orbital moments 
vary strongly with the lattice site and their signs follow the spin 
moments so that they partly cancel each other like the spin moments 
and magnetic dipole moments do. In addition, the NPs will show 
higher values of the Fe orbital moments due to surface effects. For the 
200 nm thick polycrystalline magnetite film which served as a bulk 
reference the VT was observed with VSM at 120 K as shown in the 
Supplementary Information. With field- cooled/zero -field- cooled 
measurements peaks in the blocking temperature distributions were 
found at 120 K, which indicate the VT because their temperature 
positions were independent of the magnitude of small magnetic 
fields. Interestingly, at this temperature no transition was observed 
with XMCD in accordance with the results in Ref [37]. Surprisingly, 
with XMCD a transition was detected well above the VT temper- 
ature, i.e. between 150 K and 175 K. An error in the temperature 
measurement can be excluded because with the used experimental 
setup, i.e. the high-field endstation at beamline UE46-PGM1 at the 
electron storage ring BESSY II, various transitions between 4 K and 
300 K have been found at the known transition temperatures. 
Therefore we believe that the enhanced absolute values of the mag- 
netic dipole moments of the monoclinic structure persist above the 
VT temperature until a temperature between 150 K and 175 K. This 
is in agreement with the unusual temperature dependence of the 
magnetocrystalline anisotropy that can be explained by a local charge 
and orbital ordering well above the VT as suggested by Reznicek 
et al.^^. In this sense, our work gives the first experimental evidence 
for this interpretation of changes in the electronic structure above the 
VT. 

For the NPs, the spectral shapes in Fig. 3a and Fig. 3b indicate that 
a charge transfer could be present to ligands which have not been 
washed out completely after the synthesis^^, or that the NPs are 
further oxidized at the surface^^. According to our experience both 
processes, charge transfer to ligands and oxidation at the surface, 
result in similar spectra^^. However, the transition observed with 
XMCD disppeared after storing the NPs 17 months in solution. 
The transition was re-observed with XMCD after treating the NPs 
in-situ with a Hydrogen plasma which resulted in perfect magnetite 
composition according to the XANES and XMCD spectra. Therefore 
we confidently associate the transition observed with XMCD with 
magnetite. 

Compared to the film sample, the transition temperature of 
the NPs is reduced. A reduction of the VT temperature of NPs is 
known from other experiments with different methods^^'^^ and was 
explained with spin canting and reduced thermal stability at the 
surface^^. In our NPs, spin canting has also been experimentally 
observed as a finite slope at high magnetic fields in magnetization 
hysteresis measurements with XMCD^^ and with MoEbauer spectro- 
scopy as described in Ref [45]. Thus, a reduced VT temperature 
seems to be reasonable and can explain the reduced transition tem- 
perature measured with XMCD. For the case of NPs, magnetometry 
data may not be sufficient to detect changes in the magnetic prop- 
erties caused by a phase transition, due to a strong influence of 
blocking effects as discussed in the Supplementary Information 
where we present VSM data of our NPs. This shows the importance 
of measurements of the magnetic dipole moment as a sensitive 
detector for local changes of the electronic structure. 

One might ask whether it is reasonable to attribute the observed 
transition of the magnetic dipole moment to quasiparticles called 
trimerons. Trimerons consist of a linear unit of three Fe ions on 
octrahedral sites in which a minority- spin electron from the cent- 
ral Fe^^ donor site is distributed over two adjecent Fe^^ acceptor 
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sites^^. This distribution tends to shorten the two corresponding 
Fe^^-Fe^^ distances which was recently observed with x-ray scatter- 
ing^^. Intuitively, one might expect that this donated t2g minority- 
spin electron creates an anisotropy of the spin density distribution 
which is accompanied by a magnetic dipole moment. According to 
our electronic structure calculations the B1-B3 bond length is 
shorter than other bond lenghts in the unit cell which corresponds 
to B3-B1-B3 trimerons as shown in Fig. 1. The first potential donor 
site, Bl Fe^^, is part of a trimeron and its d^z and dy^ orbitals exhibit a 
magnetic dipole moment of 0.72 fis- Contrary, the second potential 
donor site, B4 Fe^^, which carries the largest magnetic dipole 
moment of — 1.44 /Hb in its dx2 -yi orbital, is not part of a trimeron. 
Therefore it is not sufficient to assign the observed total magnetic 
dipole moment just to trimerons. 

Summary. In summary, magnetite NPs with 6 nm diameter and a 
200 nm thick polycrystalline magnetite film which served as a bulk 
reference were experimentally investigated with XANES, XMCD and 
VSM and compared with electronic structure calculations based 
on density functional theory. The measured spectra of the NPs 
offered significant changes in the white line intensity and the 
XMCD asymmetry between 50 K and 100 K. The film sample also 
showed the transition in the XMCD asymmetry which appeared 
between 150 K and 175 K. A sum rule analysis of the XMCD 
spectra revealed that the transition observed with XMCD is due to 
an increase of the effective Fe spin moment with increasing 
temperature which was attributed to a change of the dipole 
moment of the Fe 3(i spin- density distribution. This conclusion 
was verified and explained with electronic structure calculations 
based on density functional theory. The large absolute value of the 
negative magnetic dipole moment in the monoclinic structure is 
caused by the full occupancy of the d^i _yi orbital in the spin down 
channel of the Fe^^ ion on the B4 site resulting in a negative magnetic 
dipole moment, which is only partly compensated by the positive 
magnetic dipole moment of partially occupied d^^ and dy^ orbitals of 
the Fe^^ ion on the Bl site. Since the Bl site but not the B4 site is 
part of a trimeron it is not meaningful to attribute the observed total 
magnetic dipole moment just to trimerons. In the case of magnetite, 
the experimental results evidence the occurence of local charge 
and orbital ordering well above the VT as suggested in a recent 
publication^'*. In general, detailed knowledge of the magnetic 
dipole moment is a valuable tool for understanding local changes 
in the crystal- and electronic structure, not only for complex oxides 
but also complex actinide and rare-earth compounds. 

Methods 

Synthesis of nanoparticles. The NPs were prepared using a one-pot water-in-oil 
microemulsion technique as described in detail elsewhere^^. FeCl2 and FeCls were 
used as precursors, IGEPAL® CO-520 (polyoxyethylene (5) nonylphenylether) as 
stabilizing organic surfactant and ammonium hydroxide as catalyst. The water 
droplets coated by the surfactant act as nanoreactors in which the magnetite NPs are 
formed. The diameters of the NPs (6.3 ± 0.9 nm) were determined with transmission 
electron microscopy averaging over a few hundred NPs. The NPs were prepared at the 
University of Duisburg-Essen, transported in liquid solution to Berlin, deposited onto 
substrates and transferred into ultra-high vacuum for the measurements. 

X-ray absorption. The XANES and XMCD measurements were performed in the 
highfield endstation at beamline UE46-PGM1 with polarized synchrotron radiation 
from an elliptical undulator in the electron storage ring BESSY II of the Helmholtz- 
Zentrum Berlin. The total electron yield (TEY) method used in the present 
measurements to monitor the absorption of the soft x-rays has an information depth 
which is determined by the mean free path of the secondary photoelectrons of 1 nm 
to 2 nm. Since the magnetite NPs have a mean radius of 3 nm, the element- specific 
XANES and XMCD spectra represent the average over the whole particle with a 
pronounced signal from the surface with respect to the core. 

Magnetometry. The VSM measurements were performed with a Physical Properties 
Measurement System from the company Quantum Design which has a sensitivity of 
below 10"^ J/T and is operated in the laboratory cluster at the Helmholtz-Zentrum 
Berlin. 



Density functional theory. The generalized gradient approximation (GGA) as given 
by Perdew, Burke and ErnzerhoP^ plus on-site Coulomb interaction (7 to include 
strong electron correlation in the (i-orbitals of Fe were used. The Coulomb parameter 
U = 4.5 eV"^^ and the exchange parameter / = 0.89 eV"^^ were used for all ¥e-d orbitals. 
The projector augmented wave method**^'*'' as implemented in plane- wave based 
density functional code vasp^° has been used for the calculations. The plane wave 
cutoff energy was set at 400 eV energy. 12X12X12 and 12 X 12 X 4 Monkhorst- 
Packk-point grids in the Brillouin zone were used for cubic and monoclinic structures 
respectively. The geometries were optimized until the force on all each atom was 
reduced to 0.1 eV/nm. 
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